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ABSTRACT 

Dust lanes, nuclear rings, and nuclear spirals are typical gas structures in the inner region of barred galaxies. 

Their shapes and properties are linked to the physical parameters of the host galaxy. We use high-resolution 
hydrodynamical simulations to study 2D gas flows in simple barred galaxy models. The nuclear rings formed 
in our simulations can be divided into two groups: one group is nearly round and the other is highly elongated. 

We And that roundish rings may not form when the bar pattern speed is too high or the bulge central density 
is too low. We also study the periodic orbits in our galaxy models, and And that the concept of inner Lindblad 
resonance (ILR) may be generalized by the extent of X 2 orbits. All roundish nuclear rings in our simulations 
settle in the range of X 2 orbits (or ILRs). However, knowing the resonances is insufficient to pin down the exact 
location of these nuclear rings. We suggest that the backbone of round nuclear rings is the X 2 orbital family, 
i.e. round nuclear rings are allowed only in the radial range of X 2 orbits. A round nuclear ring forms exactly 
at the radius where the residual angular momentum of infalling gas balances the centrifugal force, which can 
be described by a parameter measured from the rotation curve. The gravitational torque on gas in high 
pattern speed models is larger, leading to a smaller ring size than in the low pattern speed models. Our result 
may have important implications for using nuclear rings to measure the parameters of real barred galaxies with 
2D gas kinematics. 

Subject headings: galaxies: ISM — galaxies: kinematics and dynamics — galaxies: nuclei structures — 
galaxies: hydrodynamics 


1. INTRODUCTION 


Stellar bars are the strongest disturber in the dynamical 
evolution of disk galaxies. They not only redistribute the 
stars but also have a significant impact on the interstellar 
medium (ISM) by introducing a non-axisymmetric gravi¬ 
tational torque inside the galaxy. Such a torque leads to 
the formation of interesting morphological substructures in 
the gaseous medium. Typical features include a pair of 
dust l anes related to shock s at the leading side of the bar 
(e.g., lAthanassoulal Il992bh . a circum-nuclear star-forming 
ring (or nuclear ring for brevity) near the galactic center 

(e.g..lButall986l:lGarcia-Barreto eyal.ll991l : [B^th et al.ll995l : 

iMaoz et al.l 1200 It iMazzuca et al.l l2008h . and nuclear spirals 
inside the ring which may be a cha nnel for gas infall to 
fuel the central black hole (BH) (e.g., iShlosmanet alJll990l: 
Regan & Mulchae"^ 119991 : iKnapen etak 20001; Dlaz_gL^ 


20031 Laurikainen et al.l [20041 Tvan de Ven & Fathil 120101 : 

Kim et alJ l2012alk These substructures are commonly 
observed in nearby barred g alaxies and are s tudied 


[Roberts et al.l 119791 iSchwarz 

11981: van Albada & Robertd 

1 198 It IButa & Combeslll996l 

Martini et alJ I2003albh. Thev 


may also be used as indicators to co n strain galaxy prop¬ 
erties (e.g.. IWeiner & SellwoodI 119991 IWeiner et alJ 1200 ll 
iPerez et al.ll2004l). 

Nuclear rings are thought to form as a res ult of nonlinear 
intera cti ons of gas with a ba r p otenti al ( e.g., IShlosman et aH 


19901 fildler & Shlosmanl '[19^ iPineretalJ I1995L 

Knapen et alJ 199 51: iButa & Combeslll996l I Regan & riubeni 

20031 iKim et alj|2012ai) . Due to the bar torque, gas readily 


forms dust-lane shocks at the leading side of the bar and 
flows inward along these dust lanes, forming a ring-like 


structure very close to the center. Observations show that 
usually nuclear rings are relatively circular in shape and 
have very high surface densities to trigger intense star forma¬ 
tion (e.g., iBurbidge & Burbidgel Il9^ : iKnapen et al.l 120061 
Mazzuca et al.l 120081 IComeron et alJ 2O10l: IMazzuca et alJ 

201 ll) . The typical surface densities of star formation rate 
(SFR) in nuclear rings are of order 0.1 to 1 Mq kpc“^ yr“*, 
which is about two orders of magnitude larger than the 
disk-averaged SFR densities in normal galaxies. Therefore, 
previous studies suggested that starburst activities occurring 
in nuclear rings may help to build a pseudobulge, which 
is part of the secular evolu tion of the host galaxy (e.g., 
iKormendv & Kennicuttll200^ . 

While the formation mechanism of nuclear rings is still un¬ 
der debate, it has been widely believed that the location of 
nuclear rings is directly linked to the resonance radii. For ex¬ 
ample, ICombesI (Il996h suggested that a nuclear ring forms 
near the inner Lindblad resonance (ILR) when there is only 
one ILR, as the torque outside (inside) ILR drives gas inflows 
(outflows), making it pile up at the ILR. Similarly, the nuclear 
ring forms between the inner ILR (ilLR) and outer ILR (oILR) 
when there are two ILRs. This idea of resonance-driven ring 
formation requires that the bar torque dominates the thermal 
and ram pressures of the gas both inside and outside the ILR. 
There are also some studies which li nk nuclear rings to the 
period ic or bits in the bar po tential. Sellwood & WilkinsonI 
(119931) and iSellwoodI (120141) suggested that the ring forms 
where orbits can be nested w ithout intersecting each other. 

IRegan & TeubenI (I2003L120041) argued that the ring formation 
is more deeply related to the existence of X 2 orbits rather than 
the Lindblad resonances. However, the resolution of previous 
simulations are not high enough to study the details of gaseous 
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structures in the central region of the host galaxy, and the spa¬ 
tial resolution i s very important to simulate the gas flow in 
barred galaxies dSormani et al.ll2015h . 

More recentlvT lKim et alJ d2012ah revisited the issue of the 
formation of substructures in bars using high-resolution sim¬ 
ulations by varying the gas so und speed c, and the mass of 
the central black hole (Mbh)- iKim k aO d2012al) found that 
in models with smaller Cj, nuclear rings are narrower and lo¬ 
cated farther away from the center, w hile they become m ore 
spiral-like and boarder for larger Cs- iKim et akl d2012bl) ran 
a number of numerical simulations by varying mass and ax¬ 
ial ratio of the bar to study the effect of the bar strength on 
nuclear ring formation. They found that nuclear rings form 
due to the centrifugal barrier that the inflowing gas along the 
dust lanes experiences. They showed that the ring size is in 
general smaller than the outermost ILR radius and is smaller 
for a stronger bar due to greater loss of angular momentum 
at the dust lanes, con sistent with the observational results of 
IComeron et al.l d2010h . These results imply that the location 
of the nuclear ring is not only determined by the galactic grav¬ 
itational field which generates various resonance radii and or¬ 
bital families, but also affected by the properties of the inter¬ 
stellar me dium itself such as the sound speed and ma gnetic 
field (e.g.. lKulpa-Dvbel et al.ll20 1 U lKim & Stonell2012l) . 

While the results of IKim et al.l (12012bl) are informative to 
understand why stronger bars host smaller nuclear rings, their 
models varied only the bar strength but fixed the other param¬ 
eters such as the bar patten speed fib and the axisymmetric 
part of the central galactic potential. However, Hb is one of the 
most important parameters of a bar. Hb is related to the total 
angular momentum content of the gas relative to the bar and is 
likely to control the strength and curvat ure of dust-lane shocks 
(e.g. ISanchez-Menguiano et al.|[2015h . The central galactic 
potential is thought to provide the centrifugal barrier for the 
inflowing gas. Both of these two parameters are expected to 
affect the ring siz e. In this paper, we extend the models of 
IKim et aO ( 12012bl) to explore the dependence of the ring size 
upon both Hb and the central galactic potential. To control 
the shape of the latter, we vary the bulge central density pbui 
and the Mbh- We measure the size, location, and thickness of 
the nuclear rings and study their relationships with the galac¬ 
tic parameters. These systematic experiments will offer more 
clues to the physical explanation of nuclear ring formation. 
We also compare our numerical results with observations of 
barred galaxies. 

We note a few drawbacks of our numerical models from 
the outset. First of all, we consider an infinitesimally-thin 2D 
gaseous disk embedded in a rigid stellar gravitational poten¬ 
tial, which ignores the potential effects of vertical gas motions 
introduced by finite disk thickness and back-reaction of the 
gas to the stellar components. Second, we adopt an isother¬ 
mal equation of state, which simplifies what is physically a 
much more complex, turbulent interstellar gas under radiative 
heating and cooling. Third, we neglect the gas self-gravity 
and star formation that occurs in hig h-de nsity nuclear rings. 
As we discuss the caveats in Section lATl however, these ap¬ 
proximations do not change our conclusions regarding the po¬ 
sitions and shapes of nuclear rings significantly. 

This paper is organized as follows. In Section 2, we de¬ 
scribe our galaxy models, model parameters, and the numeri¬ 
cal method. In Section 3, we present the simulation results for 
models with different fib, pbui and Mbh- We discuss the for¬ 
mation of nuclear rings and compare our simulation results 
with observations in Section 4. In Section 5, we conclude 



R [kpc] 

Figure 1. Circular velocity of gas due to each galaxy component in model 
O33C24B10 along the bar major axis (solid) and minor axis (dashed). The 
black line is the total rotation curve and the other colors represent the contri¬ 
bution of different components. The BH affects the rotation curve only in the 
very central regions (with R )S 0.1 kpc). 

with a short summary of our results. 

2. GALAXY MODELS AND NUMERICAL METHOD 

We Study how gas respond to an imposed non-axisymmetric 
bar potential using hydrodynamic simulations. The bar is as¬ 
sumed to rotate rigidly about the galaxy center with a fixed 
pattern speed fib = flbZ. We solve the dynamical equations 
in the frame co-rotating with the bar in the z = 0 plane. The 
equations of ideal hydrodynamics in this rotating frame are 




j E = -EV-u, 

(1) 

—-I-U- V 

dt J 

1 _ 2 

1 

- V<i)ext+flbR“2flb X U. 

(2) 


Here E and u denote the gas surface density and veloc¬ 
ity, respectively. Equation dTJ is the well-known continuity 
equation and Equation (|2]i is the momentum equation in co¬ 
rotating cylindrical polar coordinates. 

We use the grid - based MH D _ co de Athena 

(iGardiner & Stone l2055l: I Stone et ^ 120081: iPricel 120081 : 
IStone & Gardinei( 2009h in Cartesian geometry to integrate 
Equation ([T]i and (|2|i. Athena utilizes a higher-order Godunov 
scheme that conserves mass and momentum within machine 
precisions. It also provides several different schemes for 
integration in time, spatial reconstruction, and solution of 
the Riemann problem. We take the van Leer algorithm with 
piecewise linear reconstruction and first-order flux correction. 
Eor the Riemann solver, we use the exact nonlinear solver 
for our hydrodynamic models. Our simulation domain is 
a square box with size L = 20 kpc in each direction. We 
set up a uniform Cartesian grid with 4096 x 4096 cells 
over |Y|,|y| < L/2. The corresponding grid spacing is 
Ax = Ay = 4.8 pc. These high-resolution runs are necessary 
to explore th e pr operties of bar substructures in detail (see 
also Section 14.81) . Different from previous studies, we do 
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Figure 2. Characteristic diagram or periodic orbital families for model 
O33C24B10. The intersection of each orbit with the x-axis (the bar minor 
axis) is plotted as a function of its Jacobi constant Ej. Each calculated peri¬ 
odic orbit is given by a point and we have not joined orbits of the same family 
by a continuous curve. Red and blue dots represent stable and unstable or¬ 
bits. Note that the X 2 orbital family is stable and extends along the x-axis 
from 0.270 to 1.835 kpc. 


not impose a point symmetry relative to the galaxy center, 
thus allowing odd-m modes to grow in our models. We apply 
the outflow boundary conditions at the domain boundaries 
(i.e., at |x|, |y| = L/2). The outflow boundaries assume a zero 
gradient normal to the boundary for all flow variables except 
pressure. The solver extrapolates the required information 
from interior. 

Equation (|2]l shows that there are three parameters that af¬ 
fect the simulation results. The first is the sound speed Cs 
which effectively describes the thermal and turbulent prop¬ 
erties of gas. The second parameter is the external galactic 
potential $ext- A different shape of the potential could af¬ 
fect gas flows significantly. The third parameter is the angular 
speed of the bar fib, which is assumed to be constant over the 
radius. The setup of these three parameters is discussed in 
detail in the following. 

The real ISM is multiphase and turbulent, with temper¬ 
atures chffermgj5y_jis_jnuch_asa few o r ders of magnitude 
(e.g^ [Fiel ^ Goldsmith. & Habin^ 119691: iMcKee & Ostrikerl 
119771 ]2007lf For simplicity, we treat the gaseous disk as 
being two-dimensional, isothermal, non-self-gravitating, and 
unmagnetized. The effects of q, s elf-gravity and ma gnetic 
fields are systematic ally explored in iKim et al.l (l2012al[m and 
I Kim & Stone! (HoT^)- In this paper we fix Cs to be 8 km s“' 
and set up an initially uniform gas disk with surface density 
Eq = 10 Mq pc“^ a t all radii, which are adopted from pre - 
vious studies (e.g.. lAthanassoul^ll992bt iKirn et al.ll2012ah . 
Toomre’s Q for our gas disk is ~ 2 at all radii which means 
the disk is not strongly self-gravitating and is stable to the 
formation of spirals via gravitational instability. 

The simulations s etup for <I>ext is nearly identical to those 
in iKim et alJ (1201 2ah . Here we briefly summarize our numer¬ 
ical models and high light the difference s between the current 
models and those of iKim et akl (1201 2al) . The gaseous disk is 
placed under a $ext consisting of fo ur components: a Kuzmin- 
Toomre stellar disk (lKuzminll956h . a modified Hubble stellar 
bulge , a Ferrers ellipsoid stellar bar (iFerrersll 18871 iPfennige^ 
Il984li . and a central BH represented by a Plummer sphere. 
We emphasize that a norminal dark matter halo component 
has not been included because the disk here acts as the dark 
matte r halo to generate fla t rotation curve for historical rea¬ 
sons (lAthanassoulall 1992a) . so that we can compare our work 


directly with previous studies. These idealized models are de¬ 
signed to isolate out the effects of the pbui and Mbh on the 
formation of nuclear rings. 

To represent the elongated bar potential, we use an n = 1 
inhomogeneous Ferrers prolate spheroid: 


P = 


Pbar(l-.g7" forg<l, 

0 elsewhere, 


(3) 


g^=fla^ + (x^+z^)lb\ (4) 

The bar central density pbar. major axis a and minor axis b of 
the bar are fixed at 4.47 x 10^ Mq kpc“^, 5 kpc, and 2 kpc, 
respectively. This gives a bar axis ratio of 2.5 and a bar mass 
of 1.5 X 10'° M0. The shape and mass of the bar a re fixed for 
all of ou r models; their effects were investigated bv iKim et alJ 

(I2012bh . 

In addition to the bar, the bulge p otential also affects th e gas 
flow pattern and the ring size (e.g.. lAthanassoulall 1992bh . For 
the bulge potential, we take a modified Hubble profile; 

/ , 2 ^ 3/2 

p(0 = Pbui (^1 + ^j , (5) 


where pbui and r* repr esent the bulge central density and scale 
length, respectively. iMazzuca et al.l (120111) argued that the 
ring size is well correlated with the compactness C = Vq /Rto. 
where Rto is the turnover radius of the rotation curve that has a 
velocity Vq for the flat(ter) part. Their results show that more 
compac t (w ith smaller Rto) galaxies have smaller rings (see 
Section l43] for more details). To test this idea, it is necessary 
to consider models with different compactness values, which 
we control by varying pbui- 

Table [T] lists the model parameters and some of the simu¬ 
lations results. Column (1) gives the model names. Columns 
(2)-(4) list model parameters. Columns (6)-(9) exhibit nuclear 
ring radii along bar major and minor axes, and Column (10) 
defines the nuclear ring type (see Section |3ll. The first set 
of models in Table [T] explores the effects of Hb. while the 
models in the second set adopt a different value of pbui- In 
the third set, we vary both Hb and pbui- The last set varies 
Mbh- Figure plots the total circular rotation curves to¬ 
gether with the contribution from each component for our 
canonical model O33C24B10 with Hb = 33 km s“’ kpc”*, 
Pbui = 2.4 X 10*** Mq kpc”^, and Mbh = 1 x 10^ Mq. 

In order to avoid transients, we start with an axisym- 
metrized bar an d the bar mass is linearly ramped up to its 
final maximum (lAthanassoulall 1992bt IKirn et ^I2012al) . The 
typical time scale for bar growth here (i.e. = 2 . 1 : 

is a few hundred Myr. We have experimented with differ¬ 
ent bar growth times and confirmed that the gas distribu¬ 
tion at the end of the runs is not sensitive to the bar growth 
times, although different bar growth times may cause dis- 
ti nct initial transient gas flow p atterns, similar to the results 
in iPatsis & Athanassoulal (120001) . All simulations are run for 
5 bar revolution times, by which time the gas distribution has 
reached a quasi-steady state in the bar co-rotation frame. 


3. SIMULATION RESULTS 
3.1. Characteristic diagram of the potential 

We use characteristic diagrams to identify periodic orbital 
families in our galaxy models. The typical orbital families 
in a barred galaxy are xi, X2, and X3. The xi orbital fam¬ 
ily is the backbone of the bar, which is parallel to the bar 
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Table 1 

Model Parameters and Nuclear Ring Radii 


Model 

£2b 

Pbul 

Mbh 

Rcr/it 

Inner minor 

Outer minor 

Inner major 

Outer major 

Ring type 


km kpc“^ 

lO*" M0 kpc-3 

10*^ Mq 


(pc) 

(pc) 

(pc) 

(pc) 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(V) 

(8) 

(9) 

(10) 

O21C24B10 

21 

2.4 

10 

1.906 

715.3 

1169.4 

695.8 

979.0 

^2 

O23C24B10 

23 

2.4 

10 

1.744 

642.1 

1115.7 

676.3 

920.4 

^2 

O25C24B10 

25 

2.4 

10 

1.611 

603.0 

1032.7 

676.3 

920.4 

^2 

O27C24B10 

27 

2.4 

10 

1.500 

578.6 

1018.1 

578.6 

871.6 

^2 

O29C24B10 

29 

2.4 

10 

1.408 

520.0 

905.8 

510.2 

749.5 

^2 

O31C24B10 

31 

2.4 

10 

1.330 

485.8 

812.9 

461.4 

690.9 

^2 

O33C24B10 

33 

2.4 

10 

1.263 

437.0 

803.2 

437.0 

632.3 

^2 

O35C24B10 

35 

2.4 

10 

1.205 

393.1 

744.6 

388.2 

573.7 

^2 

O37C24B10 

37 

2.4 

10 

1.155 

349.1 

690.9 

324.7 

520.0 

^2 

O39C24B10 

39 

2.4 

10 

1.111 

329.6 

617.7 

339.3 

476.1 

^2 

O41C24B10 

41 

2.4 

10 

1.067 

349.1 

637.2 

197.7 

314.9 

^2 

O43C24B10 

43 

2.4 

10 

1.040 

173.3 

270.9 

407.7 

642.1 

Xi 

O45C24B10 

45 

2.4 

10 

1.011 

227.0 

329.6 

690.9 

1008.3 

Xi 

O47C24B10 

47 

2.4 

10 

0.984 

129.3 

197.7 

622.5 

1140.1 

ri 

O49C24B10 

49 

2.4 

10 

0.958 

148.9 

187.9 

617.6 

1066.8 

Xi 

O33C12B10 

33 

1.2 

10 

1.263 

144.0 

280.7 

549.3 

1501.4 

Xi 

O33C14B10 

33 

1.4 

10 

1.263 

119.6 

275.8 

646.9 

1652.8 

Xi 

O33C16B10 

33 

1.6 

10 

1.263 

124.5 

275.8 

588.3 

1525.9 

Xi 

O33C18B10 

33 

1.8 

10 

1.263 

280.7 

681.1 

148.9 

397.9 

X2 

033C20B10 

33 

2.0 

10 

1.263 

373.5 

734.8 

388.1 

524.9 

X2 

O33C22B10 

33 

2.2 

10 

1.263 

388.1 

788.5 

397.9 

568.8 

X2 

O33C24B10 

33 

2.4 

10 

1.263 

437.0 

803.2 

437.0 

632.3 

X2 

O33C26B10 

33 

2.6 

10 

1.263 

495.6 

837.4 

466.3 

656.7 

X2 

O33C28B10 

33 

2.8 

10 

1.263 

520.2 

891.3 

485.8 

666.5 

X2 

033C30B10 

33 

3.0 

10 

1.263 

544.4 

905.7 

520.0 

705.6 

X2 

O33C32B10 

33 

3.2 

10 

1.263 

563.9 

935.1 

524.9 

720.2 

X2 

O33C34B10 

33 

3.4 

10 

1.263 

588.4 

939.9 

539.5 

725.1 

X2 

O33C36B10 

33 

3.6 

10 

1.263 

607.9 

969.2 

544.4 

734.8 

X2 

O33C38B10 

33 

3.8 

10 

1.263 

622.5 

969.2 

573.7 

734.8 

X2 

033C40B10 

33 

4.0 

10 

1.263 

632.3 

979.0 

573.7 

739.7 

X2 

O25C16B10 

25 

1.6 

10 

1.906 

490.7 

817.8 

412.5 

681.1 

X2 

O25C32B10 

25 

3.2 

10 

1.906 

866.7 

1066.8 

700.7 

871.6 

X2 

025C40B10 

25 

4.0 

10 

1.906 

905.7 

1105.9 

754.4 

944.8 

X2 

O37C16B10 

37 

1.6 

10 

1.155 

183.1 

305.1 

568.8 

979.0 

Xi 

O37C32B10 

37 

3.2 

10 

1.155 

485.8 

832.5 

461.4 

642.1 

X2 

037C40B10 

37 

4.0 

10 

1.155 

563.9 

905.7 

529.7 

676.2 

X2 

O43C16B10 

43 

1.6 

10 

1.040 






O43C32B10 

43 

3.2 

10 

1.040 

378.4 

642.1 

378.4 

520.0 

X2 

043C40B10 

43 

4.0 

10 

1.040 

480.9 

583.5 

471.2 

583.5 

X2 

045C24B100 

45 

2.4 

100 

1.011 

21.9 

275.8 

21.9 

187.9 

X2 

045C24B400 

45 

2.4 

400 

1.011 

266.1 

549.3 

231.9 

427.2 

X2 


Note. — Simulation parameters and the resulting nuclear ring radii. Col.(l): model ID. Col.(2): pattern speed of the bar. Col.(3): bulge central 
density. Col.(4): black hole mass. Col.(5): bar co-rotation radius over bar major axis. Col.(6)-Col.(9): the inner and outer nuclear ring radii al ong 
bar major and minor axes. Col.(10): nuclear ring type. The nuclear ring radii i n the first two sets of this table is visualized in FigureslTland l 1I I 
Snapshots of gas surface density of the last two sets can be seen in Figures l^and l 121 


major axis; X 2 and are centrally concentrated and perpen¬ 
dicular to the bar major axis. Figure plots the orbital en¬ 
ergy Jacobi constant Ej vs. the distance x where each orbit 
crosses the short axis of the bar. Each sequence of periodic 
orbital families generates a continuous curve, and each dot 
represents a periodic orbit. The X 2 and x^ orbital families 
form a banana shape in this characteristic curve at small x, 
the while xi orbital family extends to a larger radius with a 
higher energy. The segment located between x = 0 and 0.3 
inside X 2 (xy) orbital family is the orbits generated by the 
BH. F or the detailed shape of thes e orbits we refer the read- 


ers to 

Combes & SandersI ( 

19811). IContoDoulos (&: GrosbOll 

1 ® 

00 

o^ 

Athanassoula (1 1992a 

). Sellwood & Wilkinson! (119931) 


andlKim et al.l (12012^. 


3.2. General evolution 

Our canonical m odel O33C24B10 is nearly identical to the 
standard model in iKim et ahl (12012ah except that q is in¬ 


creased from 5 to 8 km s“', and they used CMHOG2 code 
in cylindrical geometry. The overall evolution of both mod¬ 
els is very similar, despite the two codes being quite different, 
thus demonstrating that the simulation results are robust. 

We briefly describe the gas evolution in our canonical 
model, which can also be seen in Figure[3] At f = 0 Gyr, gas is 
distributed uniformly over the whole simulation domain. As 
the bar grows, gas begins to gather in ridges at the leading side 
of the bar. When the bar potential becomes strong enough, 
the ridges develop into sharp dust-lane shocks. These shocks 
are curved and thus prone to the wiggle instability that oc¬ 
curs due to ac cumulation of potential vorticity at the d i storte d 
shock fronts (I Wada & Kodal l2004t iKim et al]l2012aL l2014l) . 
Gas loses its angular momentum when it hits the dust-lane 
shocks, falling inward to form a nuclear ring at the position 
where the centrifugal force balances the gravity. Small dense 
clumps produced at the dust-lane shocks due to the wiggle 
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Figure 3. Time evolution of the gas surface density in the bar co-rotating 
frame in Model O33C24B10. The snapshots are in logarithm scale and the 
bar is pai'allel to y-axis. Dust lanes (shocks), nuclear rings and nuclear spirals 
are clearly visible in these snapshots. {An animation of this figure is avail¬ 
able.) 


instabi lity also fall into the ring, making it fairly inhomoge¬ 
neous (iKim et alJl2012ah . 

The bisymmetric bar potential induces weak m = 2 per¬ 
turbations in the central regions. The perturbed gas ele¬ 
ments follow slightly elliptical orbits, forming spiral struc¬ 
tures whose shape depends critically on the sign of d{n - 
Kl2)ldR, such that they are leading (trailing) where d(i2- 
K/2)/dR is positive (negative) (ICombeslll996h . Here = 
R~^d^ext/dR and = R~^ d(R^ff)/dR denote the angular 
and epicyclic frequencies (iBinnev & Tremain^ [200^. In 
our models the nuclear spiral structures appear in the early 
time, with the shape similar to the previous simulatio ns 
(e.g. lEn glmaier & ShlosmanI l2000t iMaciejewskil l2004allbt 
iThakur et al.l 2009 ). However, these nuclear spirals are de¬ 
stroyed later due to the interaction with nuclear rings. When 
t > 600 Myr, the gas flow near the nuclear ring region reaches 
a quasi-steady state, although the dust-lane shocks and nu¬ 
clear spirals decay slowly. These results are overall consis¬ 
tent with the previous low-resolution simulatio ns with a point 
symmetry relative to the galaxy center (e.g., lAthanassoulal 

ll992bl:[Englmaier & Gerh^dlll997H Regan & Teubenll2003l) . 

Evolution of other models with different Hb and pbui is qual¬ 
itatively similar except that nuclear rings are absent if Hb is 
too large or pbui is too small. Figure |4] displays snapshots of 
gas surface density in logarithmic scale at f = 800 Myr from 
the third set of models in Table [T] Note that a nuclear ring 
tends to be larger for smaller Hb and/or larger pbui- We will 
explore this trend more quantitatively in the following sec¬ 
tions. 


3.3. Quantitative analysis of nuclear rings 

Unlike nuclear spirals and dust lanes that decay over time, 
a nuclear ring is a robust feature that lasts for a long period 
of time. We run the canonical model up to 3 Gyr and confirm 
that the nuclear ring still persists, with its shape and size rela¬ 
tively unchanged. This holds true even when nuclear rings are 


subject to st ar formation and supernova feedback (ISeo & KimI 
1201 3L 120141) . demonstrating that the lifetime of nuclear rings 
could be at Gyr timescale. 

The nuclear rings in our simulations are roughly elliptical. 
Figure|5]plots the density cut profiles along the bar major and 
minor axes in the canonical model. In this paper we define 
a nuclear ring exists where the gas surface density exceeds 
100 Mq pc“^, since the gas surface density in the nuclear 
ring region of real galaxies is between 100 to 1000 Mq pc“^ 
(iKormendv & Kennicutd200^ . The surface density in the nu¬ 
clear ring region of our canonical model can be as high as 
1500 Mq pc“^. Note that the nuclear ring is denser and nar¬ 
rower along the bar major axis than the minor axis, which 
is most likely due to th e convergence of the o uter X2 orbits 
on the major axis (e.g. iRegan & TeubenI [20031) . The extent 
and thickness of the ring are shown by the arrows in Figure 
|5] We then examine the relation between the size and shape 
of the nuclear ring and the galactic parameters discussed in 
Section |2l Table [T] lists inner and outer nuclear ring radii 
along bar major and minor axes in each simulation. Note that 
model O43C16B10 does not form a nuclear ring. The ellipti¬ 
cal structure in the center of this model has a very low surface 
density, therefore it does not satisfy our nuclear ring definition 
(E > 100 Mq pc“^). Also there is no X 2 orbital family in this 
model. 

For simplicity, we name the nuclear rings according to their 
morphology. That is, the rings whose shape is similar to the X2 
orbital family are named X2-type rings', these rings are almost 
round and slightly elongated along the bar minor axis. Sev¬ 
eral examples of .r 2 -type rings can be seen in FigurelH such as 
models O33C24B10 and 043C40B10. For the nuclear rings 
highly elongated along the bar major axis which probably fol¬ 
low the xi orbital family, we name them xi-type rings (models 
O33C16B10 and O43C24B10). The definitions are the same 
as in iKim et al.l (l2012bb . 

3.4. Varying pattern speed 
3.4.1. The relation between nuclear ring size and Hb 

The pattern speed Hb is one of the most important bar pa¬ 
rameters. We first study the relation between ring properties 
and Hb in our models. We vary Hb from 21 to 49 km s“' kpc“*, 
corresponding to the co-rotation resonance radius Rcr = 1.91 a 
to 0.9643- Figure |6] plots the characteristic angular frequen¬ 
cies H-k/ 2, H, and VI + k/2 along the bar major and minor 
axes by solid and dashed lines in these models, respectively. 
The middle horizontal dotted line represents our canonical 
model with Hb = 33 km s“' kpc“* (Rcr = 1 ■2a). Strictly speak¬ 
ing, the resonance radii that we have described are only valid 
for a weak non-axisymmetric perturbation. For the strong 
bars, this estimation of the resonance radii by the linear per¬ 
turbation theory is only indicative. The bar strength parameter 
for the canonical model is 0.242 as a reference (Qi, is same for 
all models in the first set of Table [T]i, which corresponds to a 
strong ba][3 

* Wliile f2b higher than 47 1cm kpc“* con'esponds to Rcr/ a less than 
unity, which is not self-consistent due to the fact that no xj stellar orbits 
(which are the back bones of the bar) can extend beyond the co-rotation radius 
IContoDoulosll 1980h . these high ftb models are included to study the effect of 
f2b in 2 systematic way. 

^ Ql, is defined as the maximal ratio of the tangential force (mainly 
due to the non-axisymm etric bar potential) to the azimuthally-averaged ra¬ 
dial force in dfext (e .g., ICombes & Sander^ 1 198 II : ILaurikainen et alTl2004l : 
IComeron et alJ2010f) . 
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Figure 4. Gas surface density at f = 800 Myr for the models. Each panel shows the inner ±4 kpc regions where the bar is oriented vertically along the y-axis. 
Each row corresponds to models with = 1.6, 2.4, 3.2, and 4.0 X 10*® Mq kpc^^ from top to bottom, while each column is for models with fib = 25, 33, 37, 
and 43 km s“* kpc“* from left to right. The white solid lines indicate the innermost and outermost orbits. Note that model 043C16B10 does not have X 2 orbits. 


We also calculate the X 2 orbital family in these models. The 
extent of the X 2 orbital family may be used to generalize the 
concept of an ILR to finite am plitude perturbations. However, 
Ivan Albada & SandersI (Il982l) found that there are strong bars 
with no X 2 orbital family, although ILRs still exist as predicted 
by the epicyclic theory. Therefore the existence of an ILR is 
a necessary but insufficient condition for the existence of the 
X 2 orbital family. 

Figure |6] shows that changing Hb affects both ilLR and 
olLR. For Hb varing from 21 to 29 km s“* kpc“* models have 
only one ILR if we ignore the artificial cusp introduced by the 
sharp density truncation at the edge of the Ferrers bar. On the 
other hand, when Ob is larger than 31 km s“^ kpc“*, models 
have three ILRs. The distance between ilLR and oILR de¬ 


creases as Ob increases. 

Figure |2] gives the size and thickness of the rings as a func¬ 
tion of Ob in our first set of simulations. In the upper panel, 
the black solid lines draw the radii of the oILR and ilLR, while 
the green lines mark the radii of the outermost and innermost 
X 2 orbits along the bar minor axis. Note that the shape of 
nuclear rings depends on Ob. For 21 km s“' kpc~* < Ob < 
41 km s“' kpc“* (l.Ola < Rcr < L91a), the nuclear rings 
are almost circular but slightly elongated along the bar minor 
axis, so they are X 2 -type rings. The size and thickness of these 
nuclear rings along the bar minor (blue arrows) and major axis 
(red arrows) follow nearly a linear relation with Ob, such that 
a higher Ob gives a narrower and smaller ring. For these nu¬ 
clear rings, the upper ends of the blue arrows are higher than 
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Figure 5. Density cut in the canonical model O33C24B10 along bar minor 
(blue) and major (red) axes at r = 940 Myr. The dotted line indicates the 
threshold surface density E = 100 Mq pc“^ for nuclear rings. The short lines 
with arrows on both ends indicate the width of the ring in the major and minor 
axes. 



Figure 6. Frequency curves for models O21C24B10, O33C24B10, and 
O49C24B10. The solid and dashed lines represent Q —«:/2 (leftmost curves), 
Q (middle curves) and Q+k,/2 (rightmost curves) along the bar major (solid) 
and minor (dotted) axes, respectively. The horizontal dotted blue lines cor¬ 
respond to = 49,33,21 km s“^ kpc“^ from top to bottom. The vertical 
dot-dashed red line indicates the bar length a, and the upper x-axis labels 
R/a. Note that the number of ILRs depends on Qb- 

those of the red arrows, with an ellipticitjQ of eung ~ 0.1- 
0.3 at the outer boundary. The lower ends of the blue and 
red arrows are very close to each other, indicating the inner 
boundary of these rings is nearly circular. Th e observed ring 
elliptic ities are less than 0.4 as reported by iComeron et al.l 
(120101) and iMazzuca et al.l (1201 Ih . which means our numer¬ 
ical results are consistent with observations. The variation 
of nuclear ring thickness also shows a trend that the absolute 

^ Ttie ring ellipticity is defined as eri„g = 1 —hr jar, where Or and hr are the 
semi-major and semi-minor axes of the ring. 



Figure 7. Upper panel: relationships between nuclear ring size and thicicness 
with Ob. The arrows indicate the inner and outer radii of the rings along 
the bar minor (blue) and major (red) axis as defined in Figure Galaxies 
with Ob > Obcrit = 41 km s“* kpc^' have xj-type rings, while those with 
Ob < Ob,crit have X 2 -type rings. The black solid lines mai'k the positions 
of the olLR and ilLR: galaxies with Ob < 31 km s“* kpc^' have only one 
ILR. Green lines indicate innermost and outermost X 2 orbits along the bar 
minor axis. Middle and lower Panels: absolute and normalized nuclear ring 
thickness vs. Ob. The asterisks are numerical results and the solid lines are 
the fits. 

ring thickness decreases with increasing fib, but the trend is 
reversed with respect to the normalized ring thickness (the ra¬ 
tio of ring thickness to ring r adius ). We further explain and 
discuss these trends in Section|4j2] 

When fib >41 km s * kpc \ the rings suddenly become 
very elongated along the bar major axis, which can be seen 
by the large separation between the red and blue arrows (thus 
these are xi-type rings). Howev er, X] -type rings are a lmost 
never seen in observations (e.g.. IComeron et a^l2010^ . Ac¬ 
cording to this figure, we conclude that X 2 -type rings can 
only exist in our models when fib < flb.crit = 41 km s“' kpc“* 
(corresponding to Rcr = l.Qla). If fib is larger than flb.crit. 
then the rings are very elliptical along bar major axis and 
the linear correlation between fib and nuclear ring size breaks 
down. For the models with fib < flb.crit. all the nuclear rings 
are X 2 -type and their sizes correlate well with fib, although 
these models have different number of ILRs (3 ILRs for 
fib > 31 km s“^ kpc“* and 1 ILR for fib < 31 km s“’ kpc“\ see 
Figure 121). Therefore the correlation between fib and X 2 -type 
nuclear ring size is not affected by the number of the ILRs, 
which may imply that reading ILR directly off the frequency 
curves is not accurate and does not give accurate informa- 
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Figure 8. |r x F| for our Ferrers bar model. The white ellipse denotes the 
har outline. We see that the region where the largest |r X F| occurs is mainly 
at the inner edge of the bar. 


tion on the ring position, as we cautioned earlier. We remark 
that rib.crit = 41 km s“* kpc“* found here is valid only for our 
current set of model parameters; the critical patten speed de¬ 
pends on th e ba r an^alaxy parameters (for more discussions 
see Section [33] andl4li. 

Figure [T] also shows that the allowed radial range for the 
X 2 orbits follows nicely with the location of the ILRs. The 
formation of an xj -type ring is due to the limit set by the in¬ 
nermost X 2 orbit. If the bar drives gas flow inward to a radius 
that is smaller than the innermost X 2 orbit, then there are no 
stable X 2 orbits to support the inflow gas. Therefore the gas 
would instead follow the xi orbital family to remain quasi¬ 
steady, since its orbital energy is higher than that of the X 2 
orbital family for a given angular momentum and it extends 
all the way to the center (see Figure|2]l. We see that in Figure 
|4]and|2l all xi-type rings are well inside the innermost X 2 or¬ 
bit along the bar minor axis. This may explain why there is a 
break for the trend of nuclear ring sizes on Figure |7] We will 
elaborate further on this point in Section |4] 


3 . 4 . 2 . The torque exerted on the gas by the bar 

We try to explain why higher fib models tend to have 
smaller rings. The non-axisymmetric bar potential produces 
shocks in the gas flows that would otherwise follow galaxy 
rotation. Owing to a larger velocity relative to the bar poten¬ 
tial, shocks are stronger in models with relatively smaller Ob- 
Stronger shocks due to the bar potential tend to form farther 
down stream from the potential m inimum, i.e., the bar major 
axis (iHopkins & Ouat aer3l201 ll), analogous to the cases of 
spiral shocks (Kim & Ostrikeill2002L 120061: iGittins & Clarkd 
I2004t [Kim & Kimll2014l) . In Figure [8] we plot the absolute 
value of (r x F) for our canonical Ferrers bar model, where F 
is the gravitational force due to the bar. The white ellipse rep¬ 
resents the bar outline. We see that the maxima of |r x F| 
occur near the inner edge of the bar. In Figure |9] we plot 
the gas surface density E and the corresponding bar torque 


E|r X F| on gas at f = 80 Myr for models with 3 different fib 
(27,37,41 km s“' kpc“*fl The model with a higher fib pro¬ 
duces overdense ridges with a smaller offset from the bar ma¬ 
jor axis. In Figure |9] we see that dust-lane shocks in models 
with high fib tend to be located in the regions where the non- 
axisymmetric bar torque is larger, while those in models with 
low fib are located in the regions with smaller bar torque. This 
result means that the gas in higher fib models experiences a 
stronger bar torque. 

Lower fib models form larger rings because the gas loses 
less angular momentum due both to shocks and to the bar 
torque. On the other hand, models with a relatively lower 
fib have stronger shocks, but these models still apparently 
lose les s ang ular momentum in total (also see discussions in 
Section 14.21 . This suggests that the bar torque may play a 
more important role in removing angular momentum than the 
shocks. We will study the process of the angular momentum 
removal in greater depth in the future. 

3 . 5 . Varying central density 

The bulge component defines the Hubble sequence of the 
host disk galaxy, and it plays an important role in secular evo¬ 
lution. N-body simulations further sugge st that the bulge and 
bar can also interact with each other (e.g.. lAthanassoulal2002l 
IShen & Sellwoodll2004h . Therefore, the influence of bulge 
mass or bulge type/proflle on the gas flow in barred galaxies 
should also be considered. In this study we experiment with 
the impact of the bulge central density pbui on nuclear rings. 
We vary pbui from 1.2 x 10^° to 4.0 x 10^° Mq kpc“^. In or¬ 
der to keep the total mass of the bulge within 10 kpc constant 
(Mbui within 10 kpc is fixed to 3.37 x 10'° Mq, which is about 
36% of the total bulge-tbar mass within 10 kpc), the scale 
length of bulge Rb is changed accordingly from 0.467 kpc to 
0.306 kpc, with smaller Rb corresponding to a more concen¬ 
trated bulge. The variation of pbui and Rb barely affects the 
bar strength, giving Qb = 0.246 for model O33C12B10 and 
Qb = 0.234 for model 033C40B10. 

Figure[T0]shows the frequency curves in these models. Dif¬ 
ferent pbui generates a different number of resonances, sim¬ 
ilar to the case of varying Hb- In Figure [TT] we show the 
variations of the ring size and thickness with pbui- Again, 
the black solid lines indicate the location of the oILR and 
ilLR. The ring type depends not only on Hb but also on pbui- 
We And that lower pbui prefers xi-type rings and higher pbui 
tends to form larger X 2 -type rings with a critical bulge cen¬ 
tral density of pbui.cnt = 1-8 x 10'° Mq kpc“^. This result 
suggests that we need a relatively massive and concentrated 
bulge to support the formation of X 2 -type rings, which is con¬ 
sistent with previous simulatio n results (e.g., lAthanassoulal 
Il992bt I Regan & TeubenI 120031) . According to the critical 
model (O33C18B10), if the enclosed bulge mass at 1 kpc is 
less than ~ 15% of the total enclosed mass at 5 kpc, then it 
is hard to form an X 2 -type ring. The absolute thickness of nu¬ 
clear rings is not sensitive to pbuh both along the bar major and 
minor axis. In fact the absolute thickness exhibits a slightly 
negative correlation with pbuh as shown in the middle panel 
of Figure [TT] The normalized ring thickness shows a slightly 
stronger negative correction with pbui- These correlations may 
be examined with future observations. 

The correlation of pbui with the nuclear ring size becomes 
flat on the high pbui end, which may imply that pbui higher 

^ In these models the growth time of the har is fixed at 100 Myr in order to 
facilitate the comparison. 
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Figure 9. Gas surface density S (upper panels) and the corresponding bar torque E|r x F| on the gas (lower panels) for 3 different Ob models 
(27,37,41 km kpc“^). The white ellipse denotes the bar outline. We see that dust-lane shocks (overdense ridges) in models with high Qb t^nd to be lo¬ 
cated in the regions where the non-axisymmetric bar torque is large. This means that gas in higher f^b models experiences a lai'ger bar torque compai'ed to those 
in lower f^b models, leading to higher angular momentum loss in gas. 


R/a 

0.01 0.10 1.00 10.00 



R [kpc] 

Figure 10. Frequency cuiwes for models O33C12B10 (blue) and 
033C40B10 (black). Similar to Figure the solid and dashed lines rep¬ 
resent $1 —/c/2 (leftmost curves), Q, (middle curves), and Q + k,/2 (rightmost 
curves) along bar major and minor axes. When we change pbui, we also mod¬ 
ify the scale length of bulge to keep the total mass within 10 kpc constant, so 
that different pbui only affects the central region {R kpc). The horizontal 
dotted blue line con'esponds to $lb = 33 km s“^ kpc“^ for these models. Note 
that pbui changes the number of ILRs. 

than 4.0 x 10'° Mq kpc“^ generates nuclear rings similar to 
that in model 033C40B10. We conclude that .r 2 -type rings 
can only exist in our models when pbui > Pbui.cnt = 1.8 x 
10*° Mq kpc“^. Of course, pbui,crit may be different if we use 


other types of bulge prohles. Considering the effects of the 
Ob, this critical value can be smaller with a smaller Ob than 
i^b.crit = 41 km s“* kpc“* as we do see an .r 2 -type ring in model 
O25C16B10 (Figure nil, and the critical value of Ob can be 
higher with higher pbui than /Obui.cnt- Again we see that the 
correlation between pbui and X 2 -type ring size is not affected 
by the number of ilLRs either, which may be also due to the 
inaccurate estimation of resonances for a strong bar. 

The innermost and outermost X 2 orbits are plotted as green 
lines in Figure [TT] The break between the xi- and X 2 -type 
rings here is still due to the boundary set by the innermost X 2 
orbit. If we decrease pbui below pbui.crit. the X 2 orbits can no 
longer exist as the center is not dense enough. The radius of 
the innermost X 2 orbit increases with decreasing pbui. and an 
.ri-type ring may form when the gas falls inside the innermost 
X 2 orbit. 

3.6. Effects of a massive central “black hole" 

iKim et aO (12012al) has demonstrated that adding a mas¬ 
sive BH does not greatly affect the shape and size of X 2 - 
type rings, so we investigate the influence of the BH in those 
galaxies with xi-type rings. We use models O45C24B10, 
045C24B100 and 045C24B400 together to demonstrate our 
results. Figure [12] shows that model O45C24B10 forms an 
xi-type ring with outer semi-major and semi-minor axes of 
600 pc and 150 pc. The gas forms an X 2 -type ring with an 
averaged radius of ^ 200 pc in model 045C24B100 and of 
^ 500 pc in model 045C24B400. The corresponding fre¬ 
quency curves and density distributions of these three mod¬ 
els are shown in Figure [T2| We emphasize that the BHs here 
are unrealistically massive; in real galaxies the mass of the 
BH is about ~ 0.15-0.5% of the bulge mass Mbui accord- 
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Figure 11. Upper panel: relationships between nuclear ring size and thick¬ 
ness with different pbui- The an'ows indicate the inner and outer radii of the 
rings along the bar minor (blue) and major (red) axis as defined in Figure 
O Models with pbui > pbui,crit = 1-8 X 10^® Mq kpc“^ tend to have an X2- 
type ring, while those with less central density tend to have an xi -type ring. 
The black solid lines mai'k the locations of the oILR and ilLR, the second of 
which ceases to exist when pbui > 2.6 x 10^^ Mq kpc“^. Green lines indicate 
the innermost and outermost X 2 orbits along the bar minor axis. Note that the 
leftmost model in the upper panel (model O33C12B10) does not have the X 2 
orbital family. Middle and lower Panels: absolute and normalized nuclear 
ring thickness vs. pbui- The asterisks are the numerical results and the solid 
lines are the fits. 


ing to the Mnn-Mh,,! relation (e.g.. 

Masorrian et all 1 1 9981 

iGiiltekin et al.ll2009llKormendy & Hoi 

2013b. The bulge mass 


within 2 kpc here is 4.65 x 10^ Mq so that the BH should have 
a mass similar to 1-3x10^ Mq. Since we use a much more 
massive BH than 0.5%Mbui, the “BH” actually acts as another 
ultra-compact bulge in practice. 

Still, our “BH” potential dominates only the very inner part 
in the galaxy, as revealed in Figure [T] and fT2t d). The main 
difference among these models with different Mbh is the pres¬ 
ence/absence of an ilLR (or innermost .^2 orbit). By removing 
the ilLR, models 045C24B100 and 045C24B400 produce 
.r 2 -type rings, indicating that the formation of xi-type rings 
may be related to the presence of an ilLR or the existence of 
an innermost X 2 orbit, consistent with previous results. How¬ 
ever, varying the mass of the “BH” barely changes the loca¬ 
tion of oILR, as shown in Figure [T2( d). This result shows that 
one can generate quite different nuclear rings while keeping 
oILR u nchanged, which is a counter-example to the simula¬ 
tions in ISheth et ^ (|2000|) who found that the radius of the 
nuclear ring was proportional to the radius of the oILR. There¬ 
fore at least in some situations the size of nuclear rings does 


not rely only on the location of oILR. Combining with the re¬ 
sults of the last subsection, we conclude that the locations of 
2 C 2 -type rings may not be well predicted by the ILRs directly 
re ad off the frequency curves, whic h is similar to the f i nding s 
in Ivan de Ven & Chang! (1200^ and iPifiol-Ferrer et al.l (1201411 . 
Figure [ 12 ] shows that a central “BH” with proper mass could 
help to form ultra-compact nuclear rings (Rung/a < 0.1). 

We have shown previously that in order to form an X 2 -type 
ring in the central part, we need X 2 orbits to extend all the 
way to the center (in other words, we need the radius of the 
innermost X 2 orbit to be as small as possible). This can be 
achieved by a massive “BH” potential. As we see in Figure 
fT2l a massive “BH” removes the inner boundary of the X 2 or¬ 
bital family. In such a galactic potential we could form tiny 
2 i: 2 -type rings without worrying about the formation of xi-type 
rings. 

4. DISCUSSION 

4.1. The formation mechanism of nuclear rings 

The explanation for the formation of a nuclear ring is still 
under debate. It has often been argued that the location of 
a nuclear ring should be linked to the ILRs, since between 
the co -rotation radiu s and ILRs the potential drives mass in¬ 
ward (ICombeslll996h . However, at exactly which radius (rel¬ 
ative to the ILRs) a nuclear ring forms is still not well un¬ 
derstood. As Equation (O suggests, the gas flow is not only 
determined by the external gravitational held $ext and the pat¬ 
tern speed Hb, which govern all the resonances , but also af¬ 
fected by the sound speed Cj. iKim et akl (12012al) have shown 
that a larger Cj tends to yield a smaller nuclear ring, as larger 
thermal pressure bends the gas streamlines more strongly, in¬ 
ducing stronger dust-lane shocks and causing a larger amount 
of angular momentum loss. This implies a more fundamen¬ 
tal interpretation, which does not rely only on resonances, is 
needed to explain the formation of nuclear rings. 

We propose that there are three conditions to form a 
roundish X 2 -type ring. First, the galaxy needs to have enough 
gas that can lose angular momentum in the center. This may 
explain why gaseous nuclear rings are rare in early type dis k 
galaxies, as they may haye run out of gas (e.g.. lKnat)enl2005l) . 
Second, the X 2 orbital family should exist in the galaxy, as it 
is the backbone of the nuclear ring. Third, exactly where the 
ring settles is determined by the interplay of the residual an¬ 
gular momentum of the gas after interacting with the potential 
which is goyernedby gas dynamics, and the radial distribution 
of centrifugal force (or X 2 orbits) which is determined by the 
mass distribution of the galaxy. We explain the third point in 
more detail below. 

To estimate the total angular momentum remoyal of gas, 
we simplify the problem with some assumptions. We assume 
that the gas around the bar end Rbafl is nearly unperturbed at 
the beginning and on nearly circular orbits. After the impact 
of shocks and the bar torque, gas loses angular momentum 
and flows inward, producing a ring-like structure at Rung- The 
gas in the nuclear ring also follows nearly circular motions. 
iKim et aP (l2012allbh showed that a nuclear ring forms where 
the centrifugal force of the inflowing gas balances the grayity 
for models with fixed Hb and pbui- To see if this is also the 
case for models yarying Hb and pbui, we define the normalized 


^ In our models, = a and Rung is the ring radius along bar minor axis. 
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Figure 12. Gas flow patterns for (a) 045C24B10 with Mbh = 1 x 1 O’ M© , (b) 045C24B100 with Mbh = 1 x 10* M© , (c) 045C24B400 with Mbh = 4 X 10* M© 
at f = 950 Myr, and (d) the corresponding frequency curves. In (a), the white ovals plot the innermost and outermost X 2 orbits. In (b) and (c), only the outermost 
JC 2 orbit is given, since the ;C 2 orbits extend to the center due to a massive “BH”. In (d), the red, blue and black curves give — ^2/2 for models with O45C24B10, 
045C24B100, and 045C24B400, respectively. The horizontal dotted blue line marks Qb = 45 km s~^ kpc“^ for these models. Note that it is very hard to identify 
whether the nuclear structure in panel (b) and (c) is an jC 2 -type ring or a pair of tightly-wound nuclear spirals. 


centrifugal force relative to the bar end at 7? = R^ar as 

/(R) = F(J^)/Fb^ = . (6) 

I'bar/'^bar 

Here F{R) is the centrifugal force at radius R, Fbai is the 
centrifugal force at R = Rbar. and Vbai is the circular veloc¬ 
ity at R = Rbar which is read off the rotation curve. Suppose 
that the rotation curve is nearly flat with a constant veloc¬ 
ity Vc0, then the centrifugal force at the nuclear ring radius 

* The assumption of flat rotation curve is indeed simplified, but not criti¬ 
cal to our argument. Actually most ;C 2 -type roundish rings settle on the ris- 


is RCRring) = VoVf^ring OC Rjiag, while the angular momentum 
loss is ~ Vo(Rbar ^ring) — where Fq — V()Rbai> 

then R'(Riing) cx; (Lq “ AL)“' . The more angular momentum 
loss AL there is, the larger F is; therefore, F may serve 
as a rough proxy of AL. The total angular momentum re¬ 
moval during the evolution is reflected by the normalized cen- 

ing pai't of the rotation curve. Without loss of generality, we can assume 
that the velocity at R is v{R) oc R^, thus allowing one to easily show that 
F{R) oc a rotation curve generally rises from 

nearly solid-body rotation, then declines smoothly to the flat{ter) outer part. 
If a < 0.5, F{R) is a monotonically decreasing fu nctio n of R, and AL is 
still positively correlated with F as shown in Figure [Isl and the discussions 
therein. 
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R [kpc] R [kpc] 


Figure 13. Relation between nuclear ring radius and the normalized centrifugal force profile /(R). The black line in panel (a) represents the f{R) curve of the 
models with varying Q]), while the red and blue f{R) curves in panel (b) represent the models that have xi- and JC 2 -type rings with varying pbui- The open black 
squares give the numerical results for the jC 2 -type rings, while the long-dashed purple lines are their fits. The black arrow in panel (a) indicates the position of 
^peak- ^peak is generally smaller than Rto- 


trifugal force on the nuclear ring f^ng = f(Rnng) = ^ring/^bai, 
where Fnng = F{Rnng) = v{Rnng)'^/Rdng- ^ring Correlates with 
the amount of angular momentum loss AL due to shocks and 
the bar torque, and F\,^ is a normalizing constant. 

To study the distribution of the centrifugal force, we plot the 
prohle of f{R) along bar minor axis in Figure [TS] The black 
squares in panel (a) and (b) are the X 2 -type ring radii along 
bar minor axis from Figure 171 and fTTI As we assume that gas 
follows circular motions inside the nuclear ring, those squares 
should lie on the curves plotted depending on how small the 
rings are. The dashed lines ht (and extrapolate) the squares. 

We denote /?peak as the radius where f(R) is peaked, as indi¬ 
cated by the black arrow in Figure fOl' al. Since f(R) oc 
for v(R) cx /?“, /?peak Separates the steeply rising inner part of 
the rotation curve {a > 0.5) and the more slowly rising part 
in transition to the flat(ter) outer region (a < 0.5), but it is 
not identical to the turnover radius of the rotation curve, 
(usually loosely dehned as where a becomes 0). 

Note that all the roundish rings are positioned to the right 
side of /?peak, 1-6. at the decreasing part of the f{R) curves, 
indicating that a smaller X 2 -type ring loses more angular mo¬ 
mentum and has a larger f^ag- We also see that the smallest 
X 2 -type ring is close to /?peak- Since the potential does not pro¬ 
vide enough centrifugal force for the gas in circular motions 
at /? < Rpeak, one cannot generate X 2 -type rings smaller than 
our critical model (i.e., O41C24B10 and O33C18B10) sim¬ 
ply by increasing fib or decreasing pbui- This is more evident 
in Figure fT3? bl where the blue and red curves represent the 
models with X 2 - and xi-type rings, respectively. An X 2 -type 
ring cannot form when the dashed line does not intersect with 
the curve of f{R) before /?peak- When a more massive “BH” is 
added, on the other hand, it raises the inner part of f(R) and 
provides a sufficient centrifugal force barrier for an X 2 -type 
ring at radii that would otherwise form an xi-type ring (see 
Figure fT2li. If the centrifugal force is insufficient to support 
the circular motion, then the gas changes its orbits to eccen¬ 


tric ones as it still has non-zero angular momentum (therefore 
an xi -type ring). This is an alternative explanation for Figure 
11 in iKim et al.l (1201 2bh where they argue that xi-type rings 
formed due to the velocity of gas are too high to settle on the 
X 2 orbital family, and is also an alternative explanation for the 
X 2 orbit extent theory we discussed in Section [3] This expla¬ 
nation is simpler as we do not need to compute orbits. 


4.2. Qualitative explanation for the trends of nuclear ring 
size and thickness 

Here we attempt a qualitative physical explanation for the 
trends of nuclear ring size and thickness (Figure IT] and fTTT) . 
The reason that low Hb simulations form larger rings is be¬ 
cause gas loses less angular momentum compared with high 
Hb runs, as the bar torque is weaker in low Hb simulations than 
the high Hb ones. For low Hb simulations, the bar growth time 
is larger and the gas response time is longer for slow bars; the 
gas crossing the shock region at later times loses more angu¬ 
lar momentum than gas at earlier times as the shocks and bar 
torque become stronger with time. Thus the nuclear ring is 
more spatially extended in simulations with a low Hb simula¬ 
tion than high Sib. 

High pbui simulations also have sightly smaller gas response 
time since they have more mass in the central part, so higher 
Pbui simulations result in slightly narrower rings. Since the bar 
strength of these simulations is nearly the same, i.e. f^ng is al¬ 
most constant, the purple dashed line is almost horizontal as 
is seen in panel (b) of Figure [T3] High pbui simulations could 
form larger nuclear rings at a larger radius since they provide 
the “required” centrifugal force barrier hrst, as gas flows in¬ 
ward from outer part to the center and intersects hrst with the 
curves provided by higher pbui- We emphasize that the angu¬ 
lar momentum loss in different pbui simulations is nearly the 
same, making f^ng nearly constant for this set of simulations 
in which the bulge is only an axisymmetric component. 
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Figure 14. Compaiison of our numerical results with the observational mea¬ 
surements for the relation between the normalized ring radius faring^^d 
bar strength Q^. The open blac k squares represent the non-self-gravitating 
models from IKim et all j2Q12bD . the symbols give the mean values aver¬ 
aged over t = 0.3-1.0 Gyr, while error bars represe nt the standard devia¬ 
tion s. The star symbols a nd filled triangles are from IComeron et HI 420101) 
and IMazzuca et all 420111) . respectiv ely. Yellow s tars stand for the nuclear 
rings in double baiTed galaxies from lErv^ 42004) and blue stars are single 
barred or unbarred galaxies (for unbarred galaxies Qb is Qg which is defined 
in the same way). The Qj, measured in the unbarred galaxies is defined in the 
sam e way but is mainly c ontributed by the strong spiral arms or oval disks 
(see lComeron et aT]l201QI) . The red line with arrows on both ends indicates 
the range of the nuclear ring size in our models, with the largest ring gen¬ 
erated from model O21C24B10 and the smallest ring generated from model 
045C24B100. 

4.3. Relations between ILRs andx 2 orbits 

As Ivan Albada & SandersI (119821) suggested, the existence 
and extent of the X 2 orbital family can be used to general¬ 
ize the concept of an ILR for strong bars. We do see that 
the trends of the X 2 orbital family and ILRs follow each other 
closely in Figure|2]and[TT] We suggest that the innermost and 
outermost X 2 orbits represent the idea of the ilLR and oILR, 
respectively, and that they are more accurate than the linear 
perturbation theory for strong bars. Although knowing the lo¬ 
cation of resonances (or X 2 orbits) is insufficient to pin down 
the exact location of the nuclear ring, it is still informative to 
constrain the ring size. The exact location should be deter¬ 
mined by the third condition discussed above. 

4.4. An Empirical relation between /^ng and bar parameters 

The dashed lines in Figure [13] indicate the extrapolation of 
ring locations with their “required” centrifugal force barrier to 
form an X 2 -type ring. These lines that fit our numerical results 
are given by 

/(^ring) = -5.23(Rring/kpc)-l- 10.46 for varying Ob, (7) 
and 

/W-mg) = -0.04(Rring/kpc)-l-7.09 for varying pbui- (8) 

Based on the discussion above, /ring correlates with the total 
angular momentum removal due to the bar potential, which 


should be a function of Qb, a/b, and Rcs./a, but less sensi¬ 
tive to pbui as Figure [TST bl demonstrates. The reason is that 
the first three parameters influence the shock formation and 
angular momentum loss, while pbui just affects the angular 
momentum distribution in the inner part. 

For an external galaxy with a bar, a nuclear ring and a well- 
defined rotation curve, we can simply compute /ring by Equa¬ 
tion (|6]l as we know the location of the ring and where the bar 
ends. Qb or a/b can be derived from the near-lR images of 
the galaxy (e.g. IComeron et alJl2()T(ih . Qb, a/b, and Rcr/a 
all contribute to gas an gular momentum removal, which is re¬ 
flected by /ring (SectionjrTj). If we find a correlation between 
/ring and these bar parameters, we can put some constraints on 
the bar pattern speed Rcr/a for real galaxies. We explore the 
relation between /ring and Qb, ajb, and Rcr/ a in a follow-up 
paper, and apply it to a sample of observed rings (Li, Shen & 
Kim 2015, in preparation). We also suggest that since nuclear 
rings may not tightly correlate with the resonance radii read 
directly off the frequency curves, /ring would be quite useful 
to study the properties of real galaxies that contain nuclear 
rings. 

4.5. Comparison with observational results 

Observations show that almost all the nuclear rings are 
roundish (i. e. X 2 -tvpe), and they appear to be located close 
to Rio (e.g., iButa & CombesI [19961 : iPinol-Ferrer et al.l 12014 
some example galaxies are NGC1097, NGC7217, NGC5728, 
NGC1512, NGC1068, NGC3504, and ESO 565-11). Our nu¬ 
merical results are consistent with observations as most of the 
X 2 -type rings in the simulations are within 0.5 -1.0 kpc, com¬ 
parable to Rio at ~ 1.5 kpc. 

Eurthermore, Rpeak is where the rise of the inner rotation 
curve becomes shallower than since f{R) and 

Rio is where a becomes 0 as shown in Section [4T| Rpeak is 
closely related to, but smaller than RtcO- Strictly speaking, we 
argue that the site of the smallest roundish ring is constrained 
b y Rnsak instead of Rtp, a s shown in EigurelO] 

iMazzuca et al.l (1201 Ih recently measured the ring radii for 
a sample of 13 barred/unbarred galaxies that contain star¬ 
forming nuclear rings. They argued that the ring size is anti¬ 
correlated with the compactness C = Vq/Rio of the galaxy, 
where Vq is the velocity at the flat(ter) part of the rotation 
curve. Their results show that more compact (with smaller 
Rto) galaxies have a smaller ring. This result seems incon¬ 
sistent with our simulations at first glance. Eirst of all, if 
we treat the slope of the initial steeply rising component of 
the rotation curve (which is proportional to the central den¬ 
sity of the galaxy po) as the true compactness, C should 
be defined as C' = Vg/Rf^ instead, but we find that it still 
follows a negative correlation with the ring size even with 
this updated definition. The range of C in our models is 
from C - 2.96 x 10^ (km s-')^ kpc”' (model O33C12B10) to 
C- L08x 10^(kms-i)2kpc-' (model O33C40B10), compa¬ 
rable to the range covered by their sample. The discrepancy is 
proba bly because their sam ple is well-controlled. The galax¬ 
ies in IMazzuca et al.l (1201 Ih exhibit a positive correlation be¬ 
tween C and Qb, as well as an obvious negative correlation 
between the ring thickness and C, whereas Qb is nearly con¬ 
stant in our Eigure [TT] The anti-correlation between C and 
Rring in their results may simply be a reflection of a more in- 

^ It is theoretically possible that for some bulge profiles, /(r) may decrease 
monotonically, for which Itpeak cannot be defined. 
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tri nsic negative correl ation between R^ag and Qt,, as suggested 
bv iKim et alJ (1201 2bl) . 

Figure fT4l updates the comparison of the results from 
IComeron et al.l (120101) (star sym bols) with s i mulatio n results 
(squares with error-bars) from iKim et aO (I2012bl) on the 
Riingla-Qb p lane. There are two main changes in this figure 
compared to iKim et ahl (12012bl) . First, we divide the sam¬ 
ple into single barred/unbarred galaxie^ (blue stars) and dou- 
ble barred ga laxies (yellow stars) according to the catalog in 
lErwinl (l2004l) . If a secondary bar can further drive gas in th e 
nuclear ring to the central part (e.g. IShlosman et al.llT989h . 
one may then expect that the yellow stars should be systemati¬ 
cally lower than the blue stars. However, we observe no major 
difference between the blue stars and yellow stars, suggesting 
that secondary bars ma y not further drive gas flo w inward, 
similar to the results in iMacieiewski et al.l (120021) . Second, 
we plot the maximum and minimum sizes of the nuclear rings 
formed in our simulations as shown by the red line segment 
with arrows on both ends. The results imply that for a fixed 
bar strength, the size of nuclear rings can still vary by a factor 
of ~ 3, which is controlled by Hb, pbui and Mbh, as discussed 
in the present paper. 

A single bar potential is unlikely to explain ultra compact 
nuclear rings {Rung/a less than 0.1). We suggest the mag¬ 
netic held might play an important role in such compact rings 
since the magnetic held of equipartition strength with the ther¬ 
mal energy could make the ring size smaller by a factor of 
~ 2 (iKim & Stonel201^ . In additi on, a larger effecti ve sound 
speed can make the rings smaller (IKim et alJl2012ah . and an 
extremely compact center (such as the “BH” we use) can also 
generate compact rings, like in the model 045C24B100. 

4.6. Why xi-type rings are rare in Nature? 

Almost all nuclear rings that have been observed are X 2 - 
type, which may indicate that in real galaxies the bar does 
not rotate too rapidly and the central density is dense enough. 
These two conditions are not difficult to achieve. First, there 
is a natural upper limit for bar pattern speed {Res./a > 1), as 
the bar supporting orbits cannot extend beyond the co-rotation 
radius. Since bars can interact with bulges and dark matter ha¬ 
los by exchanging angular momentum, the bar pattern speed 
can decrea se with time as predicted by several N-body simula - 
tions (e.g.. lDebattista & Sellwoodl2()00l:lAthanassoulal2003l) . 
Second, the bar can also drive gas to the center, gradually re¬ 
sulting in an even denser center. The absence of xi-type rings 
may be the result of a self-regulated process in Nature. 

4.7. Limitations of this work 

In this work, we explore only the dynamical effects of Hb 
and pbui on the properties of nuclear rings. The self-gravity 
of gas does not make a significant effect on the ring posi¬ 
tion which is the main focus of the current work, since the 
gravitational force due to the self-gravity in the radial direc¬ 
tion becomes zero at the ring (or more precisely, at the radius 
near the middle of th e ring wher e the po tential reaches a min¬ 
imum). According to IKim et al.l (12012bl) . turning on gas self¬ 
gravity may make the ring about 5%- 10% larger for an ini¬ 
tially uniform gas disk. We also double-checked that includ- 

* A possibility miglit be tliat a bar liave existed in tlie past and is now 
dissolved but tlie nuclear rings still exist, or the nuclear ring in these unbarred 
galaxies might be due to strong spiral arms that could also remove angular 
momentum, so we consider unbarred galaxies together with single barred 
galaxies. 


ing self-gravity in our models tends to make the ring denser 
and clumpier, but the ring size changes little. 

There are also some other unaccounted physical mecha¬ 
nisms that may affect the angular momentum loss of the gas 
and may increase the scatter of ring sizes. A strong mag- 
netic held makes the ring rounder and smaller by about 20% 
(IKim & Stonell201^ . Sound speed, star-formation and feed¬ 
back may also play a role in the angular momentum transfer 
of gas, which could complicate the formation of nuclear rings 
further. 

Here we use a 2D model to study the gas dynamics, but the 
thin-disk approximation may not be valid at small R (of or¬ 
der of disk scale height, about a few hundred parsec). When 
the vertical degree of freedom is considered, the gas would 
exhibit motions along the vertical direction, mixing fluid ele¬ 
ments at different heights. Since the vertical velocities, which 
are at sonic levels at best, would be much smaller than the 
in-plane velocities induced by the bar, however, finite disk 
thickness is not expected to change the ring positions much. 
Therefore, our 2D models can be viewed as the first-order ap¬ 
proximation to what are more realistic 3D models. A simi¬ 
larity between 2D and 3D results can also be found by com¬ 
paring the morphologies of a ring and du st lanes in 3D SPH 
simulations (Figure 1 of IKim et al.l 1201 iT) with those in our 
current 2D simulations. Note that their 3D simulations have 
included self-gravity, gas cooling and heating, and supernova 
feedback. 

It is worth emphasizing that the aim in this study is not to 
attempt perfect realism, but rather to understand the physical 
mechanism of nuclear ring formation under some given, albeit 
imperfect, conditions of gas. 

4.8. Resolution effects 

Very recently, ISormani et al.l (I20L5I) studied gas flow pat¬ 
terns under a bar potential, using a numerical setup simi- 
lar to ours. Using a d i fferen t grid-based code developed by 
Ivan Albada & SandersI (119821) . they found that the size of nu¬ 
clear ring formed in their simulations decreases systemati¬ 
cally with increasing grid resolution, although their results 
appear not to be converged at their highest resolution (about 
0.4% of the bar size). Motivated by their results, we have 
made a convergence test by running two additional models 
with 1024 X 1024 and 2048 x 2048 cells for our hducial set of 
parameters. The corresponding grid spacing is Ajc = 19.2 pc 
and 9.6 pc, which are about 0.4% and 0.2% of the bar size. 
Figure [ffi] compares the snapshots of gas surface density at 
t=800 Myr from the additional models as well as our standard 
model with Ax = 4.8 pc. Although the ring is more widely 
distributed in a lower-resolution model, the averaged ring ra¬ 
dius is 620 pc, 593 pc, 590 pc for models with Ax = 19.2, 9.6, 
4.8 pc, respectively, indicating that the nuclear ring size is al¬ 
most the same for the last two resolutions. This confirms that 
our numerical results presented in the paper are not affected 
by resolution. 

5. CONCLUSION 

The main results of the present work are summarized as 
follows: 

1. Nuclear rings in hydrodynamical simulations - Nuclear 
ring regions are usually composed of high-density gas and 
massive young stars. In our numerical simulations as gas 
passes through the dust-lane shocks, it loses angular momen¬ 
tum, flows inward, and forms a nuclear ring where the cen¬ 
trifugal force balances the gravity. Nuclear rings are relatively 
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Figure 15. Convergence test for the resolution effects. The three panels show the gas density at r = 800 Myr for the canonical model with grid spacing of 19.2 pc, 
9.6 pc, and 4.8 pc from left to right, which corresponds to 0.4%, 0.2%, and 0.1% of the bar size. The white ellipse denotes the bar outline. 


long-lived and do not vary their shape much over time, thus 
they should be a good indicator to study gas morphology and 
galaxy properties. 

2. Effects of galactic parameters on nuclear ring size - 
In this paper we focus on studying the effects of varying bar 
pattern speed fib, bulge central density pbui. and black hole 
mass Mbh- The nuclear rings formed in our simulation can be 
divided into two groups: X 2 -type rings which are nearly round, 
and xi-type rings which are highly elongated along the bar 
major axis. We find that the size and thickness of X 2 -type rings 
are tightly correlated with the galactic properties we varied. 
Galaxies with low fib tend to form large and thick X 2 -type 
rings, while galaxies with high pbui form large but thin X 2 -type 
rings. A “BH” could help to form X 2 -type rings by removing 
the ilLR or the innermost X 2 orbits. These correlations may 
be examined in future observations. 

3. Nuclear ring formation mechanism - Observations show 
that nuclear rings seem to be preferably located near the 
turnover radius Rto of the galaxies. We suggest that the small¬ 
est X 2 -type ring forms near Rpeak where the normalized cen¬ 
trifugal force f(R) reaches its local maximum. Rpeak is related 
to, but smaller than Rto- We explain the transition of an X 2 - 
type ring to an xi-type ring is due to the distribution of the 
centrifugal force in the inner part or the presence of the in¬ 
nermost X 2 orbit (or ilLR). Although all X 2 -type nuclear rings 
in our simulations seem to be located well in the range of X 2 
orbits (or ILRs), knowing the resonance radii is insufficient to 
pin down the exact location of these nuclear rings. We suggest 
that an X 2 -type nuclear ring forms exactly at the radius where 
the residual angular momentum of infalling gas balances the 
centrifugal force, which can be described by a parameter /^ng 
computed from rotation curves, f -i^g can also be used to pre¬ 
dict the nuclear ring size since it serves as a rough proxy of 
the total angular momentum removal of gas, and should be a 
function of bar strength Qt, bar axis ratio a/b and bar pattern 
speed Rcr/ a. This will be carefully studied in a follow-up 
study. 
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